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Abstract. We present new near-infrared (JHK) bispectrum speckle-interferometry monitoring of the carbon star 
IRC+10216 obtained between 1999 and 2001 with the SAO 6 m telescope. The J-, and iT-band resolutions 
are 50 mas, 56 mas, and 73 mas, respectively. The total sequence of iT-band observations covers now 8 epochs 
from 1995 to 2001 and shows the dynamic evolution of the inner dust shell. The present observations show that 
the appearance of the dust shell has considerably changed compared to the epochs of 1995 to 1998. Four main 
components within a 0'/2 radius can be identified in the if -band images. The apparent separation of the two 
initially brightest components A and B increased from ~ 191 mas in 1995 to ~ 351 mas in 2001. Simultaneously, 
component B has been fading and almost disappeared in 2000 whereas the initially faint components C and 
D became brighter (relative to peak intensity). The changes of the images can be related to changes of the 
optical depth caused, for instance, by mass-loss variations or new dust condensation in the wind. Our recent two- 
dimensional radiative transfer model of IRC+10 216 suggests that the observed relative motion of components A 
and B is not consistent with the outflow of gas and dust at the well-known terminal wind velocity of 15 kms -1 . 
The apparent motion with a deprojected velocity of 19 kms -1 on average and of recently 27 kms -1 appears 
to be caused by a displacement of the dust density peak due to dust evaporation in the optically thicker and 
hotter environment. The present monitoring, covering more than 3 pulsation periods, shows that the structural 
variations are not related to the stellar pulsation cycle in a simple way This is consistent with the predictions of 
hydrodynamical models that enhanced dust formation takes place on a timescale of several pulsation periods. The 
timescale of the fading of component B can well be explained by the formation of new dust in the circumstellar 
envelope. 

Key words, techniques: image processing — circumstellar matter — stars: individual: IRC+10 216 — stars: 
mass-loss — stars: AGB and post-AGB — infrared: stars 



1. Introduction 

The carbon star IRC +10216 is a long-period variable star 
evolving along the Asymptotic Giant Branch (AGB). It is 
the nearest carbo n star known (distance ~ 110 — 150 pc; 
Crosas & Menten 1997 , Groenewegen et al. [1998 ) and the 
brightest 12 /im object outside the solar system (IRAS 
1986). A strong stellar wind has led to an almost complete 
obscuration of the star by dust. The mass- loss rate as mea- 
sured in CO rotati onal lines amounts to 2 — 5 x 1O _5 M0 /yr 
(Loup et al. 1993[ ). Detailed two-dimensional (2D) radia- 
tive transfer modeling shows that IRC +10216 had re- 
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cently suffered from an even higher mass-loss rate of ~ 
10" 4 M o /yr (Men'shchikov et al. |2001|) . Based on t he hig h 
mass-loss rate, long period of P = 649 d (Le Bertre 1992 ), 
and carbon-rich chemistry of the dust-shell, IRC+10 216 
is obviously in a v ery ad vanced stage of its AGB evolution 
(see, e.g., Blocker 1999| ). 

Interferometric near-infrared imaging of IRC+10216 
with angular resolutions of better than 100 mas has re- 
vealed that its dust shell is clumpy and bipolar, and is 
changing on a time scale of only ~lyr (Weigelt et al. 
1997| , Weigelt et al. |1998| [her eafter Paper I], Haniff & 
Buscher |1998[ Ost erbart et al. |2000| [hereafter Paper II], 
Tuthill et al. 2000p . In 1996, four components were identi- 
fied in the inner dust shell of IRC +10 216 within a radius 



2 



G. Weigelt et al.: Bispectrum speckle interferometry of IRC +10 216 



Table 1. Observational parameters. JD refers to the Julian date and to the photometric phase. Phase (maximum 
light) corresponds to JD=2449430 (Mar 18, 1994) as extrapolated from the light curve of Le Bertre dl992[ ) with 
P = 649 d. A c is the central wavelength and AA the FWHM bandwidth of the filters. TVt and TVr are the numbers 
of IRC +10 216 speckle interferograms and reference-star speckle interferograms, respectively. T is the exposure time 
per frame, S is the seeing (FWHM), p is the pixel size, and R is the resolution. In the last column, the reference stars 
are given. The observations of 1995-1998 were already presented in Paper II and are given here for completeness. 
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of 20 mas (Weigelt et al. [19971 Paper I, Haniff & Buscher 
1998| ) and were denoted as A, B, C, and D in order of de- 



creasing brightness (see Fig. [l]). On larger scales the enve- 
lope of IRC +10216 appears to be spherically symmetric 
(Mauron & Huggins |1999| , |2000[ ). Since most dust shells 
around AGB stars are known to be spherically symmet- 
ric on larger scales, whereas most proto-planetary neb- 
ulae (PPN) appear in axisymmetric geometry (Olofsson 
1996| ), it is likely that IRC +10216 has already entered 



the transition phase to PPN. This suggests that the break 
of symmetry already takes place at the end of the AGB 
evolution. So far, only a few AGB objects are known to 
show prominent asphericities of their dust shells in the 
near-infrared, and are therefore believed to have entered 
this transition phase at the end of their AGB life. This 
inc ludes , for instance, the carbon star CIT6 (Monnier et 
al. |000L and the oxygen-rich stars AFGL2290 (Gauger 
et al. P-999D and CIT3 (Hofmann et al. \Z00l\ ). 



Recently, we have developed a 2D radiat ive tr ansfer 
model of IRC +10 216 (Men'shchikov et al. |200l| ; here- 
after Paper III) which can explain many aspects of the 
nebula. In this model, the star is located at component B 
and is surrounded by an optically thick shell with bipolar 
cavities. The brightest southern component A is identified 
with radiation emitted and scattered in the optically thin- 
ner southern cavity of the dense shell. However, Weigelt 
et al ( |1998| ), Haniff & Buscher ( |1998| ), and Tuthill et al. 
(2000) argued that the star is located at the position of 
component A and that components B, C, and D are dust 
clouds. In the present paper, we attempt to explain the 
recent changes in the nebula in light of the 2D model of 
Paper III, and will comment on the alternative interpre- 
tations only in passing. 



High-resolution near-infrared monitoring of the com- 
ponents A, B, C, and D has already revealed that the 
dust shell of IRC +10 216 is rapidly evolving. The 6 m tele- 
scope speckle observations presented in Paper II cover five 
epochs between 1995 and 1998 and show that the sep- 
arations between the different components had steadily 
been increasing. For example, the distance between the 
initially brightest components A and B increased by 36% 
during 1995-1998. These results are in very good agree- 
ment with Keck telescope if -band observations obtained 
with the highest resolution to date at 7 epochs between 
1997 and 1999 by Tuthill et al. ( fc000|) . Such direct observa- 
tions of the dust-shell evolution offer an ideal opportunity 
to study the mass-loss process in the late stages of AGB 
evolution, revealing details of the dust formation process 
as well as the geometry and dumpiness of the stellar wind. 

This paper presents new near-infrared bispectrum 
speckle interferometry monitoring of IRC+10216 obtained 
with the SAO 6 m telescope in 1999-2001, adding 3 new 
epochs to the already available data of dust-shell evolu- 
tion. The present observations show that the appearance 
of the dust shell has considerably changed compared to 
the epochs of 1995 to 1998. This paper is organized as 
follows. In Sect. ||, if -band observations of the dust-shell 
evolution from 1999 to 2001, a comparison of J-, and 
if -band images of 1996/97 and 2001, and J-H, J-K, 
and H—K color images are presented. In Sect. ||, these 
observations are discussed on the basis of the general mor- 
phology, 2D radiative transfer models, and dust-formation 
models. Conclusions are given in Sect, m 



2. Observations 
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Fig. 1. if -band speckle reconstructions of IRC +10 216 for 8 epochs from 1995 to 2001. The total area is l"xl". All 
images are normalized to the brightest pixel and are presented with the same color table. North is up and east is to 
the left. 




Fig. 2. J-, H- and if -band speckle reconstructions of IRC +10 216 in March 2001. The total area is 1.6" x 1.6". All 
images are normalized to the brightest pixel and are presented with the same color table. North is up and east is to 
the left. The tick marks indicate the likely position of the central star. 



2.1. Data reduction 

New bispectrum speckle-interferometry observations of 
IRC +10 216 were carried out with the Russian 6 m 
telescope at the Special Astrophysical Observatory in 
September 1999, October 2000, and March 2001. The 
speckle interferograms were recorded with our HAWAII 
speckle camera in the J, and K bands. The speckle 
transfer function was derived from speckle interferograms 
of unresolved reference stars. Table [j] lists the obser- 
vational parameters of these measurements (1999-2001) 



complemented with the data of Paper II (1995-1998). 
Images of IRC +10 216 with resolutions of 50 mas (J), 
56 mas (iJ), and 73 mas (K) were reconstructed from 
the speckle interferograms using t he bi spectrum speckle- 
inter ferometry method ( Weige lt 1977] , Lohmann et al. 
19831 , Hofmann & Weigelt |l986|) . The modulus of the ob- 
ject Fourier transform (visibility function) was deter mine d 
with the speckle interferometry method (Labeyrie |l970 ). 
The H- and if -band images are diffraction-limited, the 
resolution in J is only slightly (14%) below the diffraction 
limit. Together with the data of Paper II, a time series of 
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8 if -band images, taken between 1995 and 2001, is now 
available. Additionally, two epochs are covered in J (1996, 
2001) and H (1997, 2001). Figures 0^ show the recon- 
structed if -band images of the innermost region of IRC 
+ 10 216 for all 8 epochs between 1995 and 2001, the most 
recent J, ii", and K images of March 2001, and a com- 
parison of the J, and K reconstructions from 1996/97 
and 2001. 



2.2. Evolution in the K band from 1995 to 2001 

The if -band images (Fig. |l|) show several compact com- 
ponents within a 200 mas radius. We denote these compo- 
nents as A, B, C, and D in the order of decreasing bright- 
ness in the 1996 image (cf. Paper I). As already shown 
in Paper II, the dust shell is dynamically evolving and 
changes its shape on a timescale of ~ 1 yr. For instance, 
the apparent separation of the two initially brightest com- 
ponents A and B increased from ~ 191 mas in 1995 to 
~ 351 mas in 2001. At the same time, B has faded and 
almost disappeared in 2000 whereas the initially faint com- 
ponents C and D have become brighter. In 2001, the in- 
tensity level of C increased to almost 40% of the peak 
intensity of A. The components A and C appear to have 
started merging in 2000. 

During these six years, the compact and bright core of 
component A has been changing its shape continuously. 
While it was almost spherically symmetric in 1995, it is 
getting more elongated in the following years. For exam- 
ple, the ratio of minor to major axis drops to ~ 0.77 dur- 
ing the years 1998 and 1999, with the major axis being 
aligned along the position angle —48° (measured at 35% 
intensity level) . The elongation of the inner core of A stays 
roughly constant. Additionally, the morphology of the im- 
age in 2001 gives evidence of the potential development of 
a bulge-like structure at PA=— 123° relative to the center 
of A, which may indicate the birth of a new component. 

The apparent separations of the various components 
with respect to A and B are shown in Figure |as a func- 
tion of time. The separations refer to the intensity maxima 
of the respective components which evolve along approx- 
imately constant position angles in the 1990s. Later, in 

2000 and 2001, the position angles change moderately. 

During the first 5 epochs the apparent separation 
between A and B increased by 23 mas/yr (Paper II). 
Interpreting this increase as real motion would lead to 
14 km/s within the plane of the sky (for a distance of 
130 pc). The corresponding least-squares fit of the 1995- 

2001 data gives 28 mas/yr (18 km/s). After 1999 com- 
ponent B became very faint and almost disappeared in 
2001. Therefore, its position in 2000 and 2001 is subject 
to larger uncertainties. It should be noted that the com- 
plete series from 1995 to 2001 can be best matched by a 
parabolic fit. This may suggest an apparent acceleration 
of ~ 5 mas/yr 2 . The fit is shown in Fig. |j. A similar ac- 
celeration term of ~ 3 mas/yr 2 was inferred by Tuthill et 
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Fig. 4. Apparent separation of the components A, B, 
C, and D with respect to A and B as function of time 
and photometric phase, resp. The lines give least-squares 
fits of the 1995-2001 data resulting in average veloci- 
ties of 28 mas/yr (A-B), 5 mas/yr (A-C, for 1995-1999), 
18 mas/yr (A-D), 8 mas/yr (B-C, for 1995-1999), and 9 
mas/yr (B-D). The statistical velocity error bars amount 
to ±1 mas/yr (A-C, A-D, B-D), ±2 mas/yr (A-B), and 
±3 mas/yr (B-C). The separations A-C and B-C in 2000 
and 2001 (shown in parentheses) are not well-defined due 
to the merging of A and C. Also shown is a parabolic fit of 
the apparent increase of A-B corresponding to an accel- 
eration of ~ 5 mas/yr 2 . Note, however, that the existence 
of such acceleration terms is not very certain judged from 
the present data. The error bar given by the pixel size (cf. 
Table |l|) is shown in the upper left corner of each panel. 



al. ( pOOO ) from Keck telescope K-band data for 7 epochs 
from 1997 to 1999. 

The positions of component D with respect to A are 
of particular interest since they are well defined for the 
complete data set: Whereas B considerably fades and C 
appears to start merging with the bright core, component 
D brightens and is well separated from A in all images. In 
Paper II, we derived a velocity of 21 mas/yr (13 km/s) for 
the increase of the apparent separation A-D during 1995- 
1998. With the new data we obtain 18 mas/yr (11 km/s) 
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Fig. 3. Left column: J-band speckle images of IRC +10 216 in April 1996 (top) and March 2001 (bottom) reconstructed 
with the same resolution of 149 mas. Middle column: i^-band speckle images of IRC +10 216 in January 1997 (top) 
and March 2001 (bottom) reconstructed with the same resolution of 73 mas. Right column: if-band speckle images 
of IRC +10 216 in January 1997 (top) and March 2001 (bottom) reconstructed with the same resolution of 90 mas. 
Contour lines are shown from 0.3 mag to 4.2 mag relative to the peak brightness in steps of 0.3 mag. North is up 
and east is to the left. The dashed lines indicate the biconical geometry of the cavities according to the 2D radiative 
transfer model of Paper III; the lines intersect at the position of component B. 



for 1995-2001, as shown in Fig. [I[ For the evolution of 
these components, there is no evidence of acceleration. 
The distance A-C has only slightly changed during the 6 
years. A linear regression yields 5mas/yr for 1995-1999. 
After 1999, A and C started merging and the separations 
do not follow the above relation anymore (see Fig. ||). 



In our radiative transfer modeling (Paper III), compo- 
nent A was identified with the southern cavity of a bipo- 
lar structure, whereas the location of the central star is 
at component B. In the lower panel of Fig. [I] we therefore 
present the apparent separations of the components C and 
D with respect to component B. Least-squares fits give 8 
mas/yr (5 km/s, for 1995-1999) for the separation increase 
of B-C and 9 mas/yr (6 km/s) for that of B-D. Due to the 
merging of A and C, the separation B-C increases strongly 
in 2000 and 2001. The position of B becomes ill-defined in 
the later epochs due to its dimming. 



2.3. Changes of the J HK images from 1996/97 to 2001 

Figure ^ shows the J, H, and K images of IRC +10 216 
obtained in March 2001. The resolution in H and K is 
diffraction-limited (56 and 73 mas, respectively) whereas 
in J 86% of the diffraction limit is reached (50 mas). While 
the brightness of the inner dust shell components (yellow 
regions) increases with increasing wavelength, the bright- 
ness of the faint surrounding nebula (red regions) stays 
roughly constant. In the J image, faint nebular structure 
on sub-arcsecond scales can be seen. The north-eastern 
arm (PA ~ 51° with respect to component A) appears 
to have a clumpy structure. Various components can be 
identified in this structure at different separations from 
the center. The innermost clump is located at a separation 
of 256 mas from the center of A, almost coinciding with 
the position of component D in the H and K images. A 
closer inspection reveals that this clump, in turn, resolves 
into two sub- structures of similar size. Beyond this double 
knot, two further clumps can be identified at a separa- 
tion of 378 mas and 522 mas, resp. The inner bright core 
A appears to show a multi-component structure in the J 
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band as well. It consists of two components separated by 
- 60 mas at PA - 150°. 

Figure || shows the J-, and if -band images of 
IRC +10 216 in 1996/97 compared to those in 2001. For a 
given band, the images are reconstructed with the same 
resolution (149 mas in J, 73 mas in H, and 90 mas in K). 
The dust-shell morphology and evolution observed in H 
are very similar to those in K. For instance, the appar- 
ent separation A-D increased from 163 mas in 1997 to 
245 mas in 2001, at the same time component A became 
elongated (axis ratio of 0.65 in 2001) and a bridge be- 
tween A and C has formed. In the J band, where one pre- 
dominantly observes scattered light, the images changed 
significantly as well. As in the other bands, component 
A became elongated along the same position angle. Three 
arms can be identified in north-western, north-eastern and 
south-western direction (see also Fig. ||). The appearance 
of the northern arms becomes much more prominent be- 
tween the two epochs. 



2.4. J—K and H—K color images 

To compute color images, the J, H, and K images ob- 
served in 2001 were reconstructed with the same resolu- 
tion of 115 mas. In our previous study (Paper II) we found 
that at previous epochs, when both components A and B 
were clearly visible, the peak separations in the J, H, and 
K images were the same within a few milli-arcseconds (at 
the same epochs). Therefore, we obtained almost iden- 
tical color images independently of whether the images 
were centered on component A or B. The color images 
are not very sensitive to relative shifts within some milli- 
arcseconds. 

The resulting J—K and H—K color images for the 
epoch of March 2001 are shown in Fig. ^| together with 
the intensity contour lines in the K band. The color maps 
are confined to regions which are brighter than 2% of the 
peak intensity in the J and K images {J—K) and 3% 
of the peak intensity in the H and K images (H—K), 
resp. The absolute magnitudes were calibrated with co- 
eval photometry kindly provided by B. Yudin. The inte- 
grated magnitudes in J, H, and K (aperture 10") are 5.98 
mag, 3.51 mag, and 0.52 mag, resp. In Fig. [|, white color 
refers to magnitude differences of 7.4 mag (J—K) and 3.9 
mag (H—K); the darkest grey refers to magnitude dif- 
ferences of 5.0 mag (J—K), and 2.7 mag (H—K). As in 
1997 (Paper II), component A turned out to be clearly the 
bluest part of the inner dust-shell region, as expected for 
a cavity illuminated by scattered light. Note that the star 
is obscured by an optically very thick dust shell (ry « 40, 
cf. Paper III) and therefore it must be located within the 
red area of the color maps. 



3. Discussion 



3.1. Dust and winds 

On larger scales, most AGB stars are surrounded by dust 
shells of spherical symmetry. This holds also for t h e dus t 
shell of IRC +10 216 (Mauron & Huggins |l999j |2000|) . 
However, on sub-arcsecond scales this symmetry appears 
to break down for giants in very advanced stages of their 
AGB evolution, with IRC +10 216 being their most promi- 
nent representative. The 2D radiative transfer modelling 
presented in Paper III strongly suggested that the cen- 
tral star of IRC +10 216 is surrounded by an optically 
thick dust shell with bipolar cavities of a full opening an- 
gle of 36° tilted with the southern lobe towards the ob- 
server (although there may exist different interpretations, 
see Sect. 0). The bright compact component A is not the 
direct light from the central star but is the southern cav- 
ity of this bipolar structure dominated by scattered light. 
According to this model, the carbon star is at the position 
of the fainter northern component B. 

This dust-shell model and the possible existence of ac- 
celeration within the cavities may be understood in terms 
of the beginning operation of interacting winds which later 
lead t o the shaping of planetary nebulae (Kwok |l982 ; 
Balick |l987|) . Within this scenario a fast stellar wind inter- 
acts with the fossil, slow AGB wind. The AGB mass-loss is 
thought to be aspherical and to take preferentially place in 
the equatorial plane. Accordingly, the polar regions show 
lower densities than the equatorial regions. Later, after a 
fast wind has developed, the less dense polar material is 
carved out leading to the formation of outflow cavities. 
This fast wind can be expected to turn on slowly at the 
very end of the AGB evolution resulting into the formation 
of biconical dust-shell structures on sub-arcsecond scales 
as observed for IRC +10 216. 

The high spatial and temporal resolution of the present 
observations reveals details of the mass-loss process even 
in the immediate vicinity of the dust condensation zone. 
The monitoring, covering more than 3 pulsation periods, 
shows that the structural variations are not related to the 
stellar pulsation cycle in a simple way. This is consistent 
with the predictions of hydrodynamical models that en- 
hanced dust formation takes place on a timescale of sev- 



eral pulsation periods (Fleischer et al. 1995). In Sect. 3.3 
simplified dust formation models are compared with the 
observations. 



3.2. Observed changes interpreted by 2D modeling 

The large number of epochs available now seems to make 
the interpretation of the sub-arcsecond structures and 
their evolution even more difficult and puzzling than be- 
fore. Here we suggest a physical picture of what presently 
is going on in this object based on the results of our 
previous 2D radiative transfer modeling (Paper III). The 
model was designed to derive the structure and physical 
properties of the star and its dusty envelope at a single 
moment, corresponding to the third epoch (January 23, 
1997) of our images. Our simplified attempt to handle the 
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Fig. 5. J—K (left) and H—K (right) color images (greyscale images) of the central region of IRC +10 216 in March 
2001, confined to regions which are brighter than 2% of the peak intensity in the J and K images (J—K color image) 
and 3% of the peak intensity in the H and K images (H—K color image), resp. The resolution is 115 mas. Contour 
lines refer to the if -band image and are shown from 0.2 mag to 3 mag relative to the peak brightness in steps of 0.7 mag. 
To construct the respective color images, the J, H, and K reconstructions were centered at the positions of maximum 
intensity. White regions refer to magnitude differences of 7.4 mag (J—K) and 3.9 mag (H—K); the regions shown in 
the darkest grey correspond to magnitude differences of 5.0 mag (J—K) and 2.7 mag (H—K). North is up and east is 
to the left. Dashed lines indicate the biconical geometry of the cavities according to the 2D radiative transfer model 
of Paper III. 



difficult problem of self-consistent time- dependent multi- 
dimensional radiative transfer modeling will be presented 
elsewhere (Men'shchikov et al. 2002, submitted). 

3.2.1. The star and its evolving environment 

Key point to the understanding of the structure and evo- 
lution of the sub-arcsecond environment of IRC +10 216 
is the knowledge of the position of the central star in the 
images. As we have suggested in the previous modeling 
(Paper III), the brightest compact peak A seen in all im- 
ages is not the direct light from the central star. The star 
is most likely located at the position of the fainter peak B, 
whereas component A is the radiation emitted and scat- 
tered in the optically thinner conical cavity of the optically 
thick bipolar dust shell. The even fainter components C 
and D in the H and K images were identified with smaller- 
scale deviations of the density distribution of the circum- 
stellar environment from axial symmetry (Fig. |]; Fig. 2 
in Paper III) . These inhomogeneities are less opaque than 
other, more regular regions of the compact dense shell. 

An alternative interpretation of the morphology of our 
images of IRC +10 216 would be that components A and 
B are the lobes of a bipolar nebula and that the star is 
located in the dark region between them. However, this 
interpretation disagrees with our radiative transfer mod- 
eling presented in Paper III. A reason for this is that the 
optically thick dust which would form the bipolar lobes 
cannot exist so close to the star, at half the distance be- 
tween the components. Thus, the star is probably located 
either at A or B. The fan-shaped morphology of compo- 
nent A in the J- and i7-band images, as well as other 



evidence presented in Papers II and III, strongly suggest 
that the star is located at the position of B. 

Having defined the stellar position at component B, we 
can now interpret the evolving appearance of IRC +10 216 
in terms of an increased mass-loss rate during the last ~ 
30 years. One of the main features of the entire 8-epoch 
sequence of imaging is the relative fading of the stellar 
component B (Figs. [l|, ^]). Other changes seen in the im- 
ages are the increasing distances of components A, C, D 
from the star (Fig. |j) and their greatly varying shapes. 
Our radiative transfer modeling (Paper III) has shown 
that the observed components are cavities in the dense 
opaque shell; therefore, the obse rved motion requires an 
additional discussion (Sect. |3.2.4 ). 



3.2.2. The near-IR images and model cavities 

In order to better visualize the relative location of the 
components, we have shown the structure (conical cavi- 
ties) of our model (Paper III) in the observed images by 
dashed lines in Figs. || and [|. The lines intersecting at 
the position of the central star and making an angle of 
46° between them, outline the biconical geometry of the 
cavities. 

The April 1996 J-band image in Fig. || best exhibits 
the bipolar geometry of IRC +10 216. Clearly visible at 
this wavelength are the fan shape of the brightest com- 
ponent A and even the scattered light from the opposite 
(northern) cavity on the far side of the dense shell. Faint 
direct stellar light is visible near the origin of the conical 
cavities. The star (component B) appears much brighter 
in the January 1997 i^-band image (Fig. |^), whereas the 
fan shape of the cavity (component A) is less prominent 
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but still well visible at this wavelength. The H and K im- 
ages in 1997 agree very well with our model presented in 
Paper III. The March 2001 H and K images are qualita- 
tively similar to the older images obtained in 1997, but 
distorted by the appearance of several fainter structures 
close to component D and by merging of the bright com- 
ponent A with component C. The faint direct stellar light 
has become more difficult to identify in H and it has been 
buried in the enhanced dust emission from all the other 
components in the March 2001 if -band image. 

The most recent J image of IRC +10 216 with a 3 times 
higher resolution (50 mas, Fig. ||) shows in much greater 
detail the fine structure of the envelope. The image reveals 
a rather isolated, faint peak at the position of the central 
star (marked in Fig. ||). The following evidence suggests 
that this peak is indeed the direct light from the star, not 
just a clump that happened to be there: (1) its position 
angle relative to A (PA « 20°) is the same as in all the 
images where the star (component B) is clearly visible, (2) 
the H image at the same epoch shows a (less isolated) peak 
at the same position, (3) the distance of approximately 347 
mas between the peak and the component A is consistent 
with the fit of Fig. |. 

One might think that the faint stellar peak should be 
much better visible in K band (Fig. ||), where optical 
depths are significantly lower compared to J. However, 
the greatly enhanced hot dust emission in K band may 
well make it completely invisible. In fact, from the con- 
tinuum of the model stellar atmosphere used in Paper III, 
we estimate that the stellar brightness in J and K bands 
is approximately the same. The dust model of Paper III 
would predict that the star is a factor of ~ 1.8 x 10 3 
brighter in K than in J band, whereas our new model of 
IRC+10216 (Men'shchikov et al., submitted), computed 
specifically for the latest epoch of March 2001, predicts a 
larger factor of ~ 1.3 x 10 4 . On the other hand, from the 
calibrated color images of Fig. | we know that the nebula 
becomes much brighter between J and K (and thus red- 
der), by a factor of ~ 1.4 x 10 4 . Therefore, it is unlikely 
that in K the stellar peak is visible better than in J band. 

3.2.3. Shapes of the near-IR color maps 

The high-resolution J—H, J—K, and H—K color images 
in Fig. U observed in March 2001 confirm our interpreta- 
tion that the bright components A, C, and D are, respec- 
tively, the cavity and smaller-scale inhomogeneities in the 
dense shell, not dense clumps of dust as it might appear 
from the images alone (Figs. |l|-|3|). In fact, all the bright 
components coincide with the blue areas in the color im- 
ages, the cavity A corresponding to the bluest spot. This is 
a natural consequence of lower optical depths along those 
directions from the star, with the "hot" stellar photons 
being scattered into the direction of the observer. This 
is illustrated in Fig. || by the dashed lines showing the 
geometry of our model (Paper III). The bluest, optically 
thinnest spot is located precisely inside the conical cavity. 



The star, obscured by ~ 40 mag of visual circumstellar 
extinction, is situated in the red area of the color images, 
as it was also in the previous epochs (Fig. 16 in Paper III). 

3.2.4. Moving dense layer or dust evaporation? 

The apparent motion of the components A, C, and D 
(Fig. p]) could be attributed either to the real radial ex- 
pansion of the opaque dense layer with several "holes" 
in the dense dust formation zone or to a displacement 
of the dust formation radius due to evaporation of re- 
cently formed dust by a hotter environment. Here, we 
analyze both processes and argue that the temperature- 
induced displacement of the dust formation zone is acting 
in IRC+10 216. For simplicity, we consider a single dust 
formation radius corresponding to the formation of carbon 
dust (cf. Paper III). 

One can explain the observed decreasing brightness of 
the star by assuming a monotonically increasing mass-loss 
rate and, hence, higher densities and optical depths of the 
wind in the dust formation zone. Higher mass loss and 
continuing condensation of new dust in the wind out of 
the gas phase increase the temperatures of the outflowing 
gas and dust due to backwarming. Increasing tempera- 
tures affect the location of the inner dust boundary of the 
envelope via dust sublimation, causing its displacement 
with a velocity that has nothing to do with the outflow 
motion of the envelope. 

Crucial to distinguishing between the real motion of a 
dense layer and the temperature-induced shift of the dust 
formation zone are estimates of the apparent velocities of 
the components A, C, and D relative to the star. For the 
assumed distance of 130 pc, the linear fit in Fig. Q gives for 
the brightest component A a velocity va ~ 18 kms -1 in 
the plane of sky. On the basis of our model (Paper III), one 
can derive for the component a deprojected radial velocity 
v r A ~ 19 kms -1 . 

Since the deprojected radial velocity is higher than the 
observed (terminal) wind outflow speed in IRC +10 216 of 
v ~ 15 kms -1 , it is unlikely that the observed changes re- 
flect just an expansion of a dense layer in which grains are 
forming. In fact, the standard picture of a stationary stel- 
lar wind predicts an acceleration of dust and gas within 
distances by a factor of ~ 2 -5 lar ger than the dust for- 
mation radius (Steffen et al. 1997 ). Due to the radiation 
pressure on dust grains, the wind velocity increases in this 
transition zone, approaching asymptotically the terminal 
velocity at larger distances. As our model (Paper III) as- 
sociates the observed components of IRC +10 216 with the 
dust formation zone, we expect that dust and gas are not 
yet fully accelerated there, i.e. that the radial outflow ve- 
locity v r < 15 kms -1 . 

However, the deprojected velocity v r A ~ 19 kms -1 is 
significantly larger than the expected wind velocity in the 
dust formation zone. Only if we assume an unlikely lower 



limit of 100 pc (Becklin et al. 1969), does the deprojected 
velocity v r A approach (from above) v ~ 15 kms -1 , which 
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is still too high for the dust formation zone. If, however, 
the actual distance to IRC +10 216 is larger than 130 pc, 
one would obtain v r ^ 19 kins" 1 and an even larger dis- 
crepancy. 

Moreover, there are reasons to believe that the acceler- 
ation depicted by the parabolic fit in Fig. [I] is real and that 
the actual apparent motion of A is now as fast as va ~ 26 
kms -1 . The corresponding deprojected velocity is then 
v r A ~ 28 kms -1 , much higher than the expected wind 
speed, for any realistic distance to IRC +10 216. Taken 
in context with increasing optical depths in the shell, this 
suggests that the observed motions are caused by the rapid 
dust evaporation due to backwarming and higher temper- 
atures in the dense environment formed by the increased 
mass loss. 

We believe that a reasonable interpretation of the ob- 
served changes in IRC +10 216 would be the following 
picture. During the recent period of increasing mass-loss 
which started ~ 20-30 years ago, a compact dense shell 
has formed around the star. The mass- loss rate was re- 
cently as high as M w 1O -4 M yr -1 (Paper III) and the 
innermost regions of the dust shell are expanding outward 
at a local wind velocity v ^ 10 kms -1 . Dust formation 
continues in the expanding material, thus increasing its 
optical depth and obscuring the central star. The optical 
depths in the polar regions remain significantly smaller 
than in the other regions of the dense shell, making the 
cavity (A) and the other components relatively brighter 
than at previous epochs. As the dense, increasingly op- 
tically thick dusty shell expands, steeply rising tempera- 
tures inside it (due to the backwarming from the steep- 
ening density front) inhibit further dust condensation and 
evaporate outflowing grains. In effect, these processes have 
been shifting recently the dust formation radius outward 
with an average velocity v r ~ 19 kms -1 (or as high as 
~ 30 kms -1 in 2001, if the apparent acceleration mea- 
sured in this work is real). One can predict that the star 
will remain obscured until M starts to drop back to lower 
values. In a few years from that moment, we could prob- 
ably be witnessing the star (B) reappearing whereas the 
cavities becoming relatively fainter. 



3.3. Dust formation models and the fading of B 

If the star is located at the position of component B, as is 
suggested by the two-dimensional radiative transfer model 
(see Sect. ^2 and Paper III), the fading of B might be 
caused by the formation of new dust along the line of sight 
towards the star. To investigate, whether such a scenario 
would be capable of explaining the observed time scale 
of the fading of B, we constructed a very simple gas box 
model. The idea is to follow the process of carbon grain 
formation in a gas element moving away from the star at 
a constant velocity. To characterize the changing thermo- 



dynamic conditions experienced by the gas element, we 
assume a power-law gas temperature stratification 



T g (r) 



Ro 



(i) 



and evaluate the gas density structure from mass conser- 
vation in a stationary, spherically symmetric configuration 

2 

(2) 



M 1 



Rq 



For the co-moving 
given by 

dr 

dt = — . 



gas element, the time coordinate is 



(3) 



In this gas element, we consider heteromolecular for- 
mation and growth of carbon grains which we compute 
according to the moment method derived by Gail & 
Sedlmayr (1988). The growth process includes reactions 
with the molecular species C, C2, C2H, and C2H2. In or- 
der to calculate the concentrations of the relevant carbon- 
bearing molecules, chemical equilibrium in the gas phase 
is assumed. We consider a carbon-rich element mixture 
with otherwise solar abundances. 

In Fig. |^ the resulting temporal evolution of the gas 
box is depicted for the following parameters: To = 3000 K 
at a radius Rq = 3.15 x 10 13 cm (corresponding to the 
Stefan-Boltzmann radius for Tq and a stellar luminosity 
of 15000 L ), a = 0.5 (i.e. an optically thin temperature 
distribution), M = lO -4 M yr -1 (reflecting a presently 
increased mass- loss rate of IRC +10216), a constant out- 
flow velocity of v = 15 kms -1 , and a carbon abundance 
(by number) of ec = 1.5 eo- 

In the upper panel of Fig. H the prescribed time evolu- 
tion of gas temperature (Eq. (jl|)) and density (Eq. (||)) in 
the gas box are shown. The second panel depicts the nucle- 
ation rate J*, i.e. the number of seed particles formed per 
second and hydrogen nucleus (solid line) and the number 
of dust grains per hydrogen nucleus (dash-dotted line). In 
the third panel, the degree of condensation, i.e. the frac- 
tion of condensible material that is actually condensed 
is shown (solid line). The dash-dotted line indicates the 
optical depth at 2.2 /im (measured from the star), that 
would build up in a stationary situation. To calculate the 
dust extinction coefficient, we assume the small particle 
limit of Mie theory. This is justified, since the mean par- 
ticle radius (a) reached after 15 P is 5 x 10 -2 /im, giving a 
maximum size parameter x = 27r(a)/A = 0.14 at the con- 
sidered wavelength of A = 2.2 /im. The optical properties 
of the grains are described by the complex refractive index 
of amorphous carbon tabulated in Preibisch et al. (1993). 
In the lower panel we plot the normalized intensity at 
2.2 /im that would emerge along the line of sight from the 
increasingly obscured radiation source (solid line). This 
intensity is calculated from the formal solution of the ra- 
diative transfer equation: 



Jo 
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Fig. 6. Temporal evolution of the gas box. The time coor- 
dinate is given in units of P = 649 d, the pulsation period 
of IRC +10216. Upper panel: gas temperature (Eq. (|l|), 
solid line, l.h.s. ordinate) and density (Eq. (H), dash-dotted 
line, r.h.s. ordinate); Second panel: dust nucleation rate J* 
(solid line, l.h.s. ordinate) and number of dust grains per 
hydrogen nucleus (dash-dotted line, r.h.s. ordinate); Third 
panel: degree of condensation f con d (solid line, l.h.s. ordi- 
nate) and optical depth at 2.2 /im (dash-dotted line, r.h.s. 
ordinate); lower panel: normalized emergent intensity at 
2.2 /xm (solid line) and observed intensity of component B, 
normalized to its maximum value (crosses). See text for 
more details. 



Here we assume thermal emission of the dust grains (B is 
the Kirchhoff-Planck function) and normalize the intensity 
to J = B(T ). 

For comparison, we also plotted in the lower panel of 
Fig. |^ the observed intensity of component B, normalized 
to its intensity in our 1997 speckle image (crosses). The 
intensity of component B fades from its maximum value 
to 17% of this value within about 1400 d, corresponding to 
2.2 pulsation periods of IRC +10216. This time interval is 
indicated in the lower panel of Fig. ^ by the vertical dashed 
lines. The comparison shows that dust condensation in 
front of the star in fact can reproduce the observed time 
scale of the fading of component B for realistic values of 
the parameters characterizing our simple toy model. 

In Table ^ the timescale r of the intensity drop to 
17% of the initial value is summarized for various pa- 



Table 2. Parameters of the gas box model and resulting 
fading time scales for an assumed outflow velocity of v = 
15 kms" 1 



a T po M C/O r 

K 10- 13 g/cm 3 M /yr 649 d 
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rameter combinations. We note a steep dependence of r 
on the C/O ratio, i.e. on the amount of material avail- 
able for dust formation. Consistent model calculations for 
dust driven outflows (e.g., Winters et al. 2000b, Arndt 
et al. 1997, Dominik et al. 1990) show, that the outflow 
velocity of a carbon-rich dust-driven wind also depends 
strongly on the C/O ratio. For the observed outflow veloc- 
ity of IRC +10216 of - 15 kms -1 , C/O ratios in the range 
(1.20 <> C/O <> 1.60) are required. Therefore, Table | im- 
plies, that the present day mass-loss rate of IRC +10216 in 
fact should have increased considerably above the "canon- 
ical" value of (a few) lO _5 M0/yr derived from CO ro- 
tational line observations probing the outer and there- 
fore, older parts of the circumstellar shell (e.g., Schoier & 
Olofsson 2001) or from CO infrared line profiles observed 
about one decade ago (e.g., Winters et al. 2000a). 

Once dust has formed in the circumstellar shell, its 
opacity would lead to a pronounced backwarming and 
therefore, to a steeper temperature gradient in the dust 
formation zone. Accordingly, in our simple gas box model, 
lower temperatures would be reached at the same density. 
As a result, dust growth becomes more efficient in this 
region and a reduced C/O ratio is sufficient to recover 
the observed fading time scale (see Table ||; compare the 
corresponding entries for a = 0.5 with a = 0.6 and 0.7). 

A more realistic treatment of the dust formation and 
mass- loss process in IRC +10216 would require the si- 
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multaneous solution of the time dependent hydrodynamic 
equations, taking into account the pulsation of the star, 
together with the equations describing the dust formation 
process and the radiative transfer problem. Such a con- 
sistent investigation will be presented in a forthcoming 
paper. 

An alternative interpretation of the fading of compo- 
nent B, together with its increasing separation from A, 
could be based on a scenario, where the star is located 
in the direction of component A, as argued in previous 
studies by Weigelt et al . fll998| ), Haniff & Buscher ( |1998| ), 
and Tuthill et al. ( pOOOj ) . In this case, component B would 
indicate thermal emission from a dust cloud, ejected from 
the star in a direction almost perpendicular to the line of 
sight. The temporal evolution of the brightness of B would 
then be the result of a competition between the initially 
increasing thermal emission due to ongoing dust forma- 
tion in the cloud and decreasing temperature of the cloud 
as it is moving away from the star. Preliminary results for 
this scenario, which are based on a consistent model cal- 
culation, are presented in Winters et al. (2002). Although 
the time scale of the fading of component B is reproduced 
quite well in this model, the calculated intensity ratio be- 
tween the expanding cloud and the star is smaller than 
the observed intensity ratio between components B and A 
by about a factor of 3. 

4. Summary 

We presented new near- infrared (JHK) bispectrum 
speckle-interferometry monitoring of the carbon star 
IRC+10216 obtained during 1999-2001 with the SAO 6 m 
telescope. The J-, and if -band resolutions are 50 mas, 
56 mas, and 73 mas, respectively. Together with the data 
of Paper II, the available if -band observations cover now 8 
epochs from 1995 to 2001 and show the dynamic evolution 
of the inner dust shell. Our images show very good agree- 
me nt wit h the images (1997-1999) reported by Tuthill et 
al. ( 2000 ). Four main components within a 0'/2 radius can 
be identified, which are surrounded by a fainter asymmet- 
ric nebula. The apparent separation of the two initially 
brightest components A and B increased from ~ 191 mas 
in 1995 to ~ 351 mas in 2001. At the same time, compo- 
nent B has been fading and almost disappeared in 2000 
whereas the initially faint components C and D became 
brighter. There is weak evidence for an accelerated appar- 
ent motion of component B with respect to A. This might 
be related to the beginning development of a fast polar 
wind or, as favored by the radiative transfer calculations, 
to rapid dust evaporation due to backwarming effects. The 
changes of the images can be related to changes of the op- 
tical depth caused, for example, by mass-loss variations or 
new dust condensation in the wind. The observed relative 
motion of components A and B with a deprojected veloc- 
ity of 19 kms -1 is most likely due to dust evaporation in 
the optically thicker and hotter environment. 

The present monitoring covers more than 3 pulsation 
periods and shows that the structural variations are not 



related to the stellar pulsation cycle in a simple way. 
This is consistent with the predictions of hydrodynami- 
cal models that enhanced dust formation takes place on 
a timescale of several pulsation periods. We have demon- 
strated, that formation of new dust along the line of sight 
towards the star can explain the observed fading time scale 
of component B, for reasonable values of the parameters 
involved in our simple gas box model. In particular, this 
dust formation calculation lends independent support to 
the previous finding that the present-day mass loss rate of 
IRC +10216 should be of the order of 10- 4 M o /yr. Further 
high-resolution observations will be most important for 
testing different views and models, and for better under- 
standing of the evolution of this complicated nebula. 
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